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We study the possibilities of parametric amplification of the gravitational perturbation during 
reheating in single-field inflation models. Our result shows that there is no additional growth of the 
super-horizon modes beyond the usual predictions. 
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As was realized in that parametric resonance instability occurs during reheating period when the inflaton field <f> 
oscillates. Since gravitational perturbation is coupled to the inflaton by the Einstein equation, it may also experience 
parametric resonance amplification during this stage. This issue has been studied in Refs j| and ||, and recently 
ON ', re-examined by Finelli and Brandenberger [Q. 

The gravitational potential <f> can be calculated by solving the linearized Einstein equation, however, in the case 
of the adiabatic perturbation with scales far outside the Hubble radius (the wavenumber k — > 0), it is convenient to 
. work with Bardeen parameter, 

Si , 2$ + !/- 1 $ _ 

l/j 3 1 + w 

t-H , In Eq.(|l|) the dot denotes the derivative with respective to time, H is the Hubble expansion rate and w — p/p is the 
0^ ■ ratio of the pressure to the density of the background. In the limit of k — > 0, the Bardeen parameter satisfies P|| 
ON, 

^(l + w)HC = Q. (2) 
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Q-i! During the stage of reheating, Eq. (||) becomes 



C = . (3) 



Recently Finelli and Brandenberger Q pointed out that when the inflaton field oscilates, = occurs periodically, 
so it is possible to have C 7^ 0. If it happens, the cosmological perturbation will undergo parametric amplification 
without violating causality. Specifically, they have considered a inflaton model with potential V(</>) = m 2 </> 2 /2 and 
solved it numerically. They found that £ is constant in time [Q. 

In this paper, we extend their work and consider a general single-field inflaton potential. We examine the evolution 
of £ during the reheating stage, and find no changes of £ in time. To begin with, we derive the equation of motion of 
the Bardeen parameter £ in perturbation theory, then we will take two specific models by numerical calculations to 
illustrate our general analytical result. 

II. PERTURBATION THEORY AND ANALYTIC ARGUMENT 

Working with the longitudinal gauge, the perturbed metric can be expressed in terms of the gravitational potential 

$ 

ds 2 = (1 + 2<P)dt 2 - a 2 (t)(l - 2<f>)dx i dx i , (4) 
where a(t) is the scale factor. The perturbed Einstein equation gives B, 
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$ + 3H$ + \— + 2(H + H 2 )} $ = k 2 U + H<p)5<f> 
a z 

k 2 

8(j> + 3H5<j> + (-z + V")6(f> = 4$0 - 2V' $ , 
a 1 

1 „ • 
$ + H$ = -k' 



(5) 
(6) 
(7) 

where k 2 = 8irG, V is inflaton potential, 6<f> is the perturbation to the field 0, and a prime denotes the derivative 
with respect to 0. 

Inserting (0) into (||) one can obtain the equation of motion for $, 



3> + (H - 24)$ + (-2 + 2iJ - 2if-C)$ = 



(8) 



To eliminate the singularities in the equation above when the inflaton field oscillates, one can make use of 
Sasaki-Mukhanov variable 



Q 



then Eq.(@) can be re- written as 
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The Bardeen parameter £ is related to Q by 



C=*Q 



In the expanding universe, the inflaton field satisfies the equation of motion, 

+ 3H + V' = . 
Differentiating (|T^) with respect to a = In a (note that d = H ) , we get 

H~ x {<!>)"• + 30 + (V" + 3H)H~ V = , 

where ( ) " ee j^. 

Since H = — k 2 2 /2 ||, and H = —n 2 4><j) , we have the following relation 

2H- 2 H4> - H- 2 H4> = . 
Substracting (|l4| ) from ([l3|), and simplifying it, we can obtain 



i +3H (i 
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Differentiating Q = C with respect to time t, we have 



(9) 
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Plugging Q, Q and Q above into Eq. (|l0|) , and making use of Eq.(|l|), finally we obtain an equation of motion for £ 



(18) 



Clearly, the solution of Eq.(|l8|) is that £ = when = (note that 7^ at the time when = 0). On the other 
hand, for = 0, Eq.© gives rise to £ = 0. We conclude that £ keeps unchanged during reheating. We should point 
out that the potential V in our analytical proof above is not specified, and our result applies for general single-field 
inflaton models. 
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III. NUMERICAL EXAMPLES 



To illustrate the analytical result in the last section, we take two specific models as examples by directly solving 
the perturbed Einstein equation with numerical calculations. The first model is V((f>) — Xcft^/A, the second one is 
a massive inflaton with self-coupled interaction^, V((f>) — m 2 </> 2 /2 + A</> 4 /4 . In the latter model two parameters 
are introduced. One is the inflaton mass to, another is the self-interaction coupling constant A. In our numerical 
calculations, for this model we take to -1 to be the units of the time and leave A/m 2 as a free parameter, at the same 
time, we only choose A/to 2 = 1 x IO^to^ 2 , Ito^ 2 , 1 x 10 3 to~ z 2 ( where m p i is the Planck mass) as illustrations. 

Fig.l and Fig. 2 present the evolution of Q and C for these two models, from which we can see that Q does not 
change over a period of time and £ does not change during the zero crossing of cb in the reheating stage (Note that 
Eqs.(||) and (||) hold only for the adiabatic perturbation with the wavenumber k — > 0, and we take k = for simplicity 
in the numerical calculations.). 



IV. DISCUSSION 



In this paper we have studied the evolution of perturbations in the inflationary cosmology and found no additional 
growth of gravitational fluctuations due to the oscillating inflaton field during reheating. Our result is valid for any 
single-field inflaton potential. For the multiple-field models Bassett et al. (?]] recently pointed out that there are 
possibilities of amplification of long wavelength perturbation. This important open issue deserves further study. 
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1 The model considered by Finelli and Brandenberger M corresponds to the limit of A — > 0. 
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Figure Captions 

Fig.l Evolution of Q and £ as a function of time in the model V(4>) — A0 4 /4. The initial condition is chosen as Q = — 1 
when (j) = 0.2m p ;. The solid and dashed lines represent the evolution of Q and ( respectively. Time is expressed in 

units of {rripi^fX) 

Fig. 2 Evolution of Q and ( as a function of time for a massive inflaton V{(j>) — m 2 (f> 2 /2 + A</> 4 /4. (a), (b), and (c) 
correspond to A/to 2 = 1 x 1CT 3 to„j 2 , Ito^ 2 , 1 x 10 3 TO pJ 2 respectively. Time is in units of to -1 . 
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Fig. 2(a) 
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Fig. 2(b) 




Fig. 2(c) 



